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We investigate whether quantum theory can be understood as the continuum limit of a mechanical theory, in 
which there is a huge, but finite, number of classical “worlds,” and quantum effects arise solely from a 
universal interaction between these worlds, without reference to any wave function. Here, a “world” means an 
entire universe with well-defined properties, determined by the classical configuration of its particles and 
fields. In our approach, each world evolves deterministically, probabilities arise due to ignorance as to which 
world a given observer occupies, and we argue that in the limit of infinitely many worlds the wave function can 
be recovered (as a secondary object) from the motion of these worlds. We introduce a simple model of such a 
“many interacting worlds” approach and show that it can reproduce some generic quantum phenomena—such 
as Ehrenfest’s theorem, wave packet spreading, barrier tunneling, and zero-point energy—as a direct 
consequence of mutual repulsion between worlds. Finally, we perform numerical simulations using our 
approach. We demonstrate, first, that it can be used to calculate quantum ground states, and second, that it is 
capable of reproducing, at least qualitatively, the double-slit interference phenomenon. 
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I. INTRODUCTION 

The role of the wave function differs markedly in various 
formulations of quantum mechanics. For example, in the 
Copenhagen interpretation, it is a necessary tool for 
calculating statistical correlations between a priori classical 
preparation and registration devices [1]; in the de Broglie- 
Bohm (dBB) interpretation, it acts as a pilot wave that 
guides the world’s classical configuration [2]; in the many- 
worlds (MW) interpretation, it describes an ever-branching 
tree of noninteracting quasiclassical worlds [3]; and in 
spontaneous collapse models, its objective “collapse” 
creates a single quasiclassical world [4]. 

In other formulations, the wave function does not even 
play a primary role. For example, in Madelung’s quantum 
hydrodynamics [5], Nelson’s stochastic dynamics [6], and 
Hall and Reginatto’s exact uncertainty approach [7], the 
fundamental equations of motion are formulated in terms of 
a configuration probability density P and a momentum 
potential S (or the gradient of the latter), with a purely 
formal translation to a wave function description via 
vp := p l /~ exp[iS/h]. These approaches can describe the 
evolution of any scalar wave function on configuration 
space, which includes any fixed number of spinless 
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particles, plus bosonic fields. In this paper, we similarly 
treat spinless and bosonic degrees of freedom. 

More recently, it has been observed by Holland [8] and 
by Poirier and co-workers [9-11] that the evolution of 
such quantum systems can be formulated without reference 
to even a momentum potential S. Instead, nonlinear Euler- 
Lagrange equations are used to define trajectories of a 
continuum of fluid elements, in an essentially hydrody- 
namical picture. The trajectories are labeled by a continu¬ 
ous parameter, such as the initial position of each element, 
and the equations involve partial derivatives of up to fourth 
order with respect to this parameter. 

In the Holland-Poirier hydrodynamical approach, the 
wave function plays no dynamical role. However, it may be 
recovered, in a nontrivial manner, by integrating the 
trajectories up to any given time [8]. This has proved a 
useful tool for making efficient and accurate numerical 
calculations in quantum chemistry [9,10]. Schiff and 
Poirier [11], while “drawing no definite conclusions,” 
interpret their formulation as a “kind of ‘many worlds’ 
theory,” albeit they have a continuum of trajectories (i.e., 
flow lines), not a discrete set of worlds. 

Here, we take a different but related approach, with the 
aim of avoiding the ontological difficulty of a continuum 
of worlds. In particular, we explore the possibility of 
replacing the continuum of fluid elements in the 
Holland-Poirier approach by a huge but finite number of 
interacting “worlds.” Each world is classical in the sense 
of having determinate properties that are functions of 
its configuration. In the absence of the interaction with 
other worlds, each world evolves according to classical 
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Newtonian physics. All quantum effects arise from, and 
only from, the interaction between worlds. We therefore 
call this the many interacting worlds (MIW) approach to 
quantum mechanics. A broadly similar idea has been 
independently suggested by Sebens [12], although without 
any explicit model being given. 

The MIW approach can only become equivalent to 
standard quantum dynamics in the continuum limit, where 
the number of worlds becomes uncountably infinite. 
However, we show that even in the case of just two 
interacting worlds it is a useful toy model for modeling 
and explaining quantum phenomena, such as the spreading 
of wave packets and tunneling through a potential barrier. 
Regarded as a fundamental physical theory in its own right, 
the MIW approach may also lead to new predictions arising 
from the restriction to a finite number of worlds. Finally, it 
provides a natural discretization of the Holland-Poirier 
approach, which may be useful for numerical purposes. 
Before considering how its dynamics might be mathemati¬ 
cally formulated and used as a numerical tool, however, we 
give a brief discussion of how its ontology may appeal to 
those who favor realist interpretations. 

A. Comparative ontology of the MIW approach 

At the current stage, the MIW approach is not yet well 
enough developed to be considered on equal grounds with 
other long-established realistic approaches to quantum 
mechanics, such as the de Broglie-Bohm and many-worlds 
interpretations. Nevertheless, we think it is of interest to 
compare its ontology with those of these better known 
approaches. 

In the MIW approach there is no wave function, only a 
very large number of classical-like worlds with definite 
configurations that evolve deterministically. Probabilities 
arise only because observers are ignorant of which world 
they actually occupy, and so assign an equal weighting to all 
worlds compatible with the macroscopic state of affairs they 
perceive. In a typical quantum experiment, where the 
outcome is indeterminate in orthodox quantum mechanics, 
the final configurations of the worlds in the MIW approach 
can be grouped into different classes based on macroscopic 
properties corresponding to the different possible outcomes. 
The orthodox quantum probabilities will then be approx¬ 
imately proportional to the number of worlds in each class. 

In contrast, the dBB interpretation postulates a single 
classical-like world, deterministically guided by a physical 
universal wave function. This world—a single point of 
configuration space—does not exert any backreaction on 
the guiding wave, which has no source but which occupies 
the entire configuration space. This makes it challenging to 
give an ontology for the wave function in parts of 
configuration space so remote from the “real” configuration 
that it will never affect its trajectory (a nice analogy can be 
found in Feynman’s criticism of classical electromagnetism 
[13]). Furthermore, this wave function also determines a 


probability density for the initial world configuration 
[2,14]. From a Bayesian perspective, this dual role is not 
easy to reconcile [15]. 

In the Everett or MW interpretation, the worlds are 
orthogonal components of a universal wave function [3]. 
The particular decomposition at any time, and the identity 
of worlds through time, is argued to be defined (at least well 
enough for practical purposes) by the quantum dynamics 
which generates essentially independent evolution of these 
quasiclassical worlds into the future (a phenomenon 
called effective decoherence). The inherent fuzziness of 
Everettian worlds is in contrast to the corresponding 
concepts in the MIW approach of a well-defined group 
of deterministically evolving configurations. In the MW 
interpretation, it is meaningless to ask exactly how many 
worlds there are at a given time, or exactly when a 
branching event into subcomponents occurs, leading to 
criticisms that there is no precise ontology [16]. Another 
difficult issue is that worlds are not equally “real” in the 
MW interpretation, but are “weighted” by the modulus 
squared of the corresponding superposition coefficients. As 
noted above, in the MIW approach, all worlds are equally 
weighted, so that Laplace’s theory of probability is suffi¬ 
cient to account for our experience and expectations. 

The skeptical reader may wonder whether it is appro¬ 
priate to call the entities in our MIW theory worlds at all. 
After all, each world corresponds to a set of positions of 
particles in real (3D) space. How, then, is the nature and 
interaction of these worlds any different from those of 
different gas species, say A and B, where the positions of all 
the A molecules constitute one world and those of the B 
molecules (each one partnered, nominally and uniquely, 
with one of the A molecules) constitute another world? The 
answer lies in the details of the interaction. 

In the above example, any given A molecule will interact 
with any B molecule whenever they are close together in 
3D space. Thus, a hypothetical being in the “A world,” 
made of some subset of the A molecules, would experience 
the presence of B molecules in much the same way that it 
would feel the presence of other A molecules. By contrast, 
as will be shown later in this paper, the force between 
worlds in our MIW approach is non-negligible only when 
the two worlds are close in configuration space. It would be 
as if the A gas and B gas were completely oblivious to each 
other unless every single A molecule were close to its B 
partner. Such an interaction is quite unlike anything in 
classical physics, and it is clear that our hypothetical 
A-composed observer would have no experience of the 
B world in its everyday observations, but by careful 
experiment might detect a subtle and nonlocal action on 
the A molecules of its world. Such action, though involving 
very many, rather than just two, worlds, is what we propose 
could lie behind the subtle and nonlocal character of 
quantum mechanics. 
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B. Outline of the paper 

In Sec. II, we introduce the MIW approach in detail. The 
approach allows flexibility in specifying the precise 
dynamics of the interacting worlds. Hence, we consider 
a whole class of MIW models which, as we discuss in 
Sec. IIB, should agree with the predictions of orthodox 
quantum mechanics in the limit of infinitely many worlds. 
We present a particularly simple model in Sec. Ill and 
demonstrate in Sec. IV that it is capable of reproducing 
some well-known quantum phenomena: Ehrenfest’s theo¬ 
rem, wave packet spreading, barrier tunneling, and zero- 
point energy. We conclude with a numerical analysis of 
oscillator ground states and double-slit interference phe¬ 
nomena in Secs. V and VI, respectively. 

II. FORMULATION OF THE MANY 
INTERACTING WORLDS APPROACH 


p,{ q) = l^(q)| 2 . (2) 


The velocity of the world-particle x is then defined by 


m k x k (t ) = 


dS r { q) 


dq k 


q=x(f) 


( 3 ) 


where 


Sf(q) = ftarg[tf ( (q)]. (4) 

This equation of motion guarantees that the probability 
density for the world configuration x ( t ) at any time t is given 
by P t (x). This property is known as equivariance [14,18]. 

In Bohm’s original formulation [2], the law of motion 
[Eq. (3)] is expressed equivalently by the second-order 
equation 


A. From dBB to MIW 

For pedagogical reasons, we introduce MIW from the 
perspective of the dBB interpretation of quantum mechan¬ 
ics, and, hence, start with a brief review of the latter. We 
regard the MIW approach as fundamental, but to show that 
it should (in an appropriate limit) reproduce the predictions 
of orthodox quantum mechanics, it is convenient to build 
up to it via the equations of dBB mechanics. It is well 
known that the dBB interpretation reproduces all the 
predictions of orthodox quantum mechanics, so far as 
the latter is well defined [14,17]. 

Consider then a universe comprising J scalar nonrela- 
tivistic distinguishable particles, in a D-dimensional space. 
(As mentioned in the Introduction, we could also include 
bosonic fields, but for simplicity we omit them here.) In 
dBB mechanics, there is a universal (or world) wave 
function \D f (q) defined on configuration space. Here, the 
argument of the wave function is a vector q = 
{q l ,..., q K } J of length K = DJ. We order these variables 
so that for D = 3, the vector (q 3k ~ 2 , r/ 3/_l , r/ 3/ ) T describes 
the position of the y'th particle. For convenience, we 
associate a mass with each direction in configuration space, 
so that the mass of the yth particle appears 3 times, as 
m 3 - 7-2 = /u 3/ ~ 1 = /« 3/ . Then, we can write Schrodinger’s 
equation for the world wave function \P f (q) as 


ih 


d 

dt 


*r(q) 


r k 


E 


2 m k 



+ V(q) 


*r(q)- 


(1) 


As well as the world wave function, there is another part of 
the dBB ontology, corresponding to the real positions of all 
the particles. We call this the world-particle, with position 
x(f) = {x'(f), ...,x K (f)} T , which we also call the world 
configuration. Reflecting “our ignorance of the initial 
conditions” [2], the initial position x(0) of the world- 
particle is a random variable distributed according to 
probability density P 0 (x), where 


m k x k = f k (x) + r k (x), (5) 

with Eq. (3) applied as a constraint on the velocity at the 
initial time (t = 0). Here, the force has been split into 
classical (f) and quantum (r) contributions, the latter called 
r because of its locally repulsive nature, which will be 
shown later. These are defined by 

f(q) = —W(q), r,(q) = -V0,( q) (6) 

(with the Arth component of V being d/dq k ). Here, 



( 7 ) 


was called the quantum potential by Bohm [2], and 
vanishes for h = 0. 

It can be shown that Eq. (5) reproduces Eq. (3) at all 
times. Although Eq. (6) looks like Newtonian mechanics, 
there is no conserved energy for the world-particle alone, 
because Q, is time dependent in general. Moreover, the 
wave function \Eh(q) evolves in complete indifference to 
the world-particle, so there is no transfer of energy there. 
Thus, these dynamics are quite unlike those familiar from 
classical mechanics. 

Suppose now that instead of only one world-particle, as 
in the dBB interpretation, there were a huge number N of 
world-particles coexisting, with positions (world configu¬ 
rations) x 1; ...,x„, ...,x N . If each of the N initial world 
configurations is chosen at random from P 0 (q), as 
described above, then 


J\)(q)»tf ‘^^((q-x^o)) (8) 

n= 1 

by construction. The approximation is in the statistical 
sense that the averages of any sufficiently smooth function 
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q>( q), calculated via either side, will be approximately 
equal, and it becomes arbitrarily good in the limit N -*■ oo. 
Clearly, a similar approximation also holds if the empirical 
density on the right-hand side is replaced by a suitably 
smoothed version thereof. 

By equivariance, the quality of either approximation is 
then conserved for all times. 

One can thus approximate P r (q), and its derivatives, 
from a suitably smoothed version of the empirical density 
at time t. From this smoothed density, one may also obtain a 
corresponding approximation of the Bohmian force 
[Eq. (6)], 

rf(q) » *jv(q; X f ) for V » l , (9) 

in terms of the list of world configurations 

X r = {xi(f),x 2 (f), ...xjv(f)} (10) 

at time t. Note, in fact, that since only local properties of 
P t { q) are required for r r (q), the approximation r, v (q; X,) 
requires only worlds from the set X, which are in the 
A'-dimensional neighborhood of q. That is, the approxi¬ 
mate force is local in /^dimensional configuration space. 

We now take the crucial step of replacing the Bohmian 
force [Eq. (6)], which acts on each world-particle x„(t) via 
Eq. (5), by the approximation r N (x„(t) m , X,). Thus, the 
evolution of the world configuration x„(f) is directly 
determined by the other configurations in X,. This makes 
the wave function \b f (q), and the functions P,( q) and S,(q) 
derived from it, superfluous. What is left is a mechanical 
theory, referred to as MIW, which describes the motion of a 
“multiverse” of N coexisting worlds Xj(t), ...x n (t),... 
x N (t), where each world configuration x n (t) is a K vector 
specifying the position of J = K/D particles. 

While the MIW approach was motivated above as an 
approximation to the dBB interpretation of quantum 
mechanics, we have the opposite in mind. We regard 
MIW as the fundamental theory, from which, under certain 
conditions, dBB can be recovered as an effective theory 
provided N is sufficiently large; see Sec. IIB. Note that the 
MIW approach is conceptually and mathematically very 
different from dBB. Its fundamental dynamics are 
described by the system of N x J x D second-order differ¬ 
ential equations 

m k x k n {t) = f k (x n (t)) + 4(x„(f);X f ). (11) 

In the absence of an interworld interaction, corresponding 
to the classical limit r^(x„(f); X,) — 0 in Eq. (11), the 
worlds evolve independently under purely Newtonian 
dynamics. Hence, all quantumlike effects arise from the 
existence of this interaction. It is shown in Secs. Ill-VI that 
this nonclassical interaction corresponds to a repulsive 
force between worlds having close configurations, leading 


to a simple and intuitive picture for many typical quantum 
phenomena. 

Note also that we use the term “MIW approach’' rather 
than “MIW interpretation.” This is because, while some 
predictions of quantum mechanics, such as the Ehrenfest 
theorem and rate of wave packet spreading, will be seen to 
hold precisely for any number of worlds, other predictions 
can be accurately recovered only under certain conditions 
in the limit N -*■ oo. This has two immediate implications: 
the possibility of experimental predictions different from 
standard quantum mechanics, due to the finiteness of N, 
and the possibility of using the MIW approach for approxi¬ 
mating the dynamics of standard quantum systems in a 
controlled manner. In the latter case, Eq. (11) must be 
supplemented by suitable initial conditions, corresponding 
to choosing the initial world configurations randomly from 
P 0 , and the initial world-particle velocities from S Q via 
Eq. (3) (see Sec. IIC below and Sec. VI). 

B. Probabilities and the quantum limit 

While each world evolves deterministically under 
Eq. (11), which of the N worlds we are actually living 
in, compatible with the perceived macroscopic state of 
affairs, is unknown. Hence, assertions about the configu¬ 
ration of the J particles in our world naturally become 
probabilistic. In particular, for a given function q>(x) of the 
world configuration, only an equally weighted population 
mean 

(<K x ))x = ^X^( x «)’ ( 12 ) 

n= 1 

over all the worlds compatible with observed macroscopic 
properties, can be predicted at any time. 

We now show that under certain conditions the MIW 
expectation values [Eq. (12)] are expected to converge to 
the ones predicted by quantum theory when the number of 
worlds N tends to infinity. For this, suppose we are 
provided a solution to the Schrodinger equation 'I', on 
A'-dimensional configuration space, and the initial con¬ 
figurations of the N worlds, x^O), x 2 (0),..., x A ,(0) are 
approximately distributed according to the distribution 
P 0 (q) — 'l'o(q ) | 2 . Hence, at t = 0 one has, for any smooth 
function cp on configuration space, 

(fl»>* 0 - J rf ql^o(q)IV(q) 

= ^( X ))xo’ 

11 = 1 

with the approximation becoming arbitrarily good as 
N —> oo for q> sufficiently regular. Suppose further that 
X V, is such that at initial time t = 0 the velocities of the 
worlds fulfill Eq. (3). As a consequence of Eq. (9), the 
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trajectory of the rath world generated by Eq. (11) should 
stay close to the corresponding Bohmian trajectory gen¬ 
erated by Eq. (5). Hence, by equivariance of the latter 
trajectories, the world configurations x 1 (t), ...,x N (t) at 
time t will be approximately distributed according to the 
distribution P t (x) = |T',(x)| 2 . That is, 

(fl»(*))*, = J ^q|'^ f (q)l 2 ^(q) 

= ( 13 ) 

n= 1 

as desired. 

In summary, the configuration-space expectation values 
of the MIW approach should coincide with those for the 
quantum state T', as N —► oo provided that at some time r 
the following conditions are met: 

(1) the worlds x^t),..., x N (r) are |4r r | 2 -distributed, 

(2) the velocities of the worlds fulfill Eq. (3) at t — r. 
Since configuration space expectation values are all that is 
required to establish empirical equivalence with orthodox 
quantum mechanics (just as in dBB [2]), this establishes the 
viability of the MIW approach. 

It is a remarkable feature of the MIW approach that only 
a simple equal weighting of worlds, reflecting ignorance of 
which world an observer occupies, appears to be sufficient 
to reconstruct quantum statistics in a suitable limit. 
Similarly, in this limit, the observer should see statistics 
as predicted by quantum mechanics when carrying out a 
sequence of experiments in his or her single world. In 
particular, from the typicality analysis by Durr et al. for the 
dBB interpretation [18], it follows that Born’s rule holds for 
dBB trajectories belonging to typical initial configurations, 
where in dBB typical stands for almost all with respect to 
the 'I', | 2 measure. Hence, since the MIW trajectories are 
expected to converge to dBB trajectories and the world 
configurations are |T',| 2 distributed in the limit described 
above, Born’s rule will hold for typical MIW worlds, where 
we emphasize again that in our MIW approach typicality 
simply means for the great majority of worlds, since each 
world is equally weighted. 

C. Which initial data give rise to quantum behavior? 

We argue above that, given a solution to the Schrodinger 
equation, one can generate corresponding initial data for 
the MIW equations of motion [Eq. (11)] whose solution 
approximates quantum theory as N -* oo. Suitable initial 
data are as per conditions (1) and (2) above. This suggests a 
converse question: given a solution, X(f) = {x 1; ...,x K } to 
the MIW equations of motion [Eq. (11)], is there is any 
solution 4', to the Schrodinger equation that can be 
approximately generated by X(?)7 

The sense in which an approximation 4', « 4', could be 
generated is the following. Given the world configurations 


X(f), one can construct approximations 4',| 2 for the 
probability density as per Eq. (8). Further, from the 
velocities X(f), approximations arg 4>, for the phase can 
be constructed via Eqs. (3) and (4). Thus, there are many 
ways a quantum mechanical candidate for 4' ( could be 
constructed. The relevant measure of the quality of the 
approximation at time t is then given by the L 2 norm 

\\%-U(t-z)%\\ 2 , (14) 

where U(t) denotes the corresponding Schrodinger evolu¬ 
tion and r is some initial time. 

From the above discussion, it might seem obvious that to 
obtain approximate quantum evolution (in this sense) one 
must simply impose the velocity constraint [Eqs. (3) and 
(4)] at the intial time t — r, for some 4' r . However, on 
further reflection, one realizes that this is no constraint at 
all. For any finite number N of worlds, there is always some 
complex function 4>(q) such that setting 4' r (q) = 4>(q) 
will match the initial velocities at the positions of those 
worlds. The point is that this 4> may fail to yield 
an approximate solution at later times, so that [| 4',— 
U(t - t) 4>|| 2 is not small, for any approximate 
reconstruction 4', and t — t sufficiently large. Thus, a more 
relevant constraint may be that the velocity constraint 
[Eqs. (3) and (4)] holds for some 4', that is smoothly 
varying on the scale of the maximum interworld distance in 
configuration space. We emphasize, however, that it is not 
completely obvious that such a constraint is necessary— 
i.e., quantum behavior may be typical in the MIW approach 
as N -* oo—a point to which we return in Sec. VII (see 
also Sec. IV). 

D. Interworld interaction potential 

The general MIW approach described in Sec. IIA is 
complete only when the form of the force r /v (x; X) 
between worlds, in Eq. (11), is specified. There are 
different possible ways of doing so, each leading to a 
different version of the approach. However, it is natural to 
seek a formulation of the MIW approach in which the 
interaction force between worlds is guaranteed to be 
conservative, that is, a force of the form 

r w (x„;X) = -V Xn U N (X), (15) 

for some potential function U N (X) defined on the N world 
configurations X = (x lt ...,x N ), where V x denotes the 
gradient vector with respect to x„. 

If we have an interworld interaction potential U N (X), 
this immediately allows the equations of motion [Eq. (11)] 
to be rewritten in the equivalent Hamiltonian form 

*« = V p ff„(X, P), p„ = -V Xn H N (X, P), (16) 

with all its attendant advantages. Here, P = (pj, ...,p„) 
defines the momenta of the worlds, with components 
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p k =jn k x k , (17) 

and the Hamiltonian is given by 

"•(X' P )=tt|? + E V '( , -) + [ '»(X). (18) 

n= 1 k= 1 n= 1 

We refer to U N (X) as the interworld interaction potential. 

To motivate the form of suitable interworld potentials 
Un, note that, to reproduce quantum mechanics in the limit 
N —► oo, it is natural to require that the average energy per 
world, (E) n = N~ l H n (X, P), approaches the quantum 
average energy in this limit. Note that we cannot associate 
a definite energy with each world, because of the interworld 
interaction. Note also that we have substituted the subscript 
N for the subscript X (or X, P as would be needed in this 
case), for ease of notation. 

Now, if the configurations x { ,...,x N sample the dis¬ 
tribution 'I'(x) | 2 , the quantum average energy can be 
written, using 'I' = P 1 ? 2 exp[iS/h], as [2] 


(£>* = / dqP{q) 


dS 


|E 2 1* ^) +i/(q) 



Here, the left-hand side is to be understood as an approxi¬ 
mation of the right-hand side, obtained via a suitable 
smoothing of the empirical density in Eq. (8), analogous 
to the approximation of the quantum force r/y(q) by 
r N (x\X) in Eq. (9). It is important to note that a good 
approximation of the quantum force (which is essential to 
obtain quantum mechanics in the large N limit), via 
Eqs. (15) and (19), is not guaranteed by a good approxi¬ 
mation in Eq. (20). This is easy to see in the case that the 
empirical density is smoothed using only members of a 
(large) fixed subset of the N worlds, implying that the force 
in Eq. (15) will vanish for all worlds not contained in this 
subset. For example, in the K — 1 case considered below, 
one could choose the subset of odd-numbered worlds (i.e., 
n odd), in which case the even-numbered worlds would 
behave classically. Thus, Eq. (20) is a guide only, and 
Eq. (9) must also be checked. 

The interworld potential JJ N in Eq. (19) is positive by 
construction. This leads directly to the existence of a 
minimum energy and a corresponding stationary configu¬ 
ration for the N worlds, corresponding to the quantum 
ground-state energy and (real) ground-state wave function 
P~ 1/2 (q) in the limit N -> oo. Note that iyf v (q; X) in 
Eq. (20) approximates Nelson’s osmotic momentum [6], 
suggesting that U N in Eq. (19) may be regarded as a sum of 
nonclassical kinetic energies, as is explored in Sec. IV. The 
quantity g k N { q; X) also approximates the imaginary part of 
the “complex momentum” that appears in complexified 
Bohmian mechanics [19,20]. 

Further connections between U N (X), g Ar (q;X), and 
r iV (x; X) are discussed in Appendix A. 

UI. SIMPLE EXAMPLE 


(for N sufficiently large). Moreover, the average energy per 
world is given via Eq. (18) as 


(E) n = N~ 1 H n (X, P) 



Assuming that the trajectories generated by the 
Hamiltonian H N are close to the dBB trajectories for 
sufficiently large N, then p k — m k x k « dS(x n )/dx k , and 
comparing the two averages shows that U N (X) should be 
chosen to be of the form 

t/ iV (X) = ^^-i T [4(x„;X)] 2 , (19) 

n= 1 k= 1 Zm 


A full specification of the MIW dynamics requires 
generation of a suitable force function r^x^X) in 
Eq. (11) or a suitable function g^( q;X) in Eq. (19). As 
discussed in the previous section, these may be viewed as 
corresponding to approximations of — VQ and VP/P, 
respectively. One may obtain many candidates for inter¬ 
world forces and potentials in this way, and it is a matter 
of future interest to determine what may be the most 
natural ones. 

In this section, we demonstrate in the simplest case of 
one dimension and one particle per world (K = D = 
7=1) how the proposed MIW approach can be math¬ 
ematically substantiated, using a Hamiltonian formulation. 
In Sec. IV, we show that this example provides a nice 
toy model for successfully describing various quantum 
phenomena. 


where 


g k N ( q;x) 


h 1 dP( q) 
2P(q) dq k 


A. One-dimensional toy model 

For N one-dimensional worlds, it is convenient to 
distinguish them from the general case by denoting the 
configuration x n and momentum p„ of each world by x n 
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and p n , respectively. It is also convenient to label the 
configurations such that x l < x 2 < ■ ■ ■ < x N (it will be 
seen that this ordering is preserved by the repulsive nature 
of the interaction between worlds). Similarly, we use X = 
(xx ,..., x N ) in place of X, and P in place of P. Note that we 
can thus consider X and P as vectors. The mass of each 
one-dimensional world-particle is denoted by m. 

The aim here is to obtain a simple form for the interworld 
potential U N in Eq. (19). It is easiest to first approximate 
the empirical distribution of worlds by a smooth probability 
density P x (q) and use this to obtain a suitable form for 
g k N {q\X) in Eq. (20). 

Now, any smooth interpolation P x (q) of the empirical 
density of the N worlds must satisfy 

1 , r 

t j^2<P(x n )~ / dqP x {q)(p{q) 

V n=\ 4 

N fx„ 

dqP x {x n )(p{x n ) 

11= 1 J X n -1 

N 

= ~ X„-l)P X (x„)<p(x n ) 

n =1 

for sufficiently slowly varying functions cp(x) and suffi¬ 
ciently large N. This suggests the following ansatz, 


Px(x n ) = 


N(x n -x n _i) N(x n+l -x n )' 


( 21 ) 


for the smoothed distribution P x at x = x n , which relies on 
an assumption that the interworld separation is slowly 
varying. The success of this ansatz requires that the 
dynamics we derive based upon it preserve this slow 
variation, which appears to be the case from the simulations 
in Sec. VI. Under this assumption, the two expressions in 
Eq. (21) have a relative difference 0(l/N). Further, for 
large N, we can approximate the derivative of P x (q) at 
q = x n by 


To ensure that the summation is well defined for n = 1 and 
n = N, we formally define x 0 = — oo and x N+ \ = oo. 

B. Basic properties of the toy model 

The MIW toy model corresponding to Eq. (23) is defined 
by the equations of motion [Eq. (16)] and the Hamiltonian 
[Eq. (18)]. Since the total energy is conserved, two adjacent 
worlds x n and x n+l cannot approach each other arbitrarily 
closely, as the potential U N would become unbounded. 
This implies a mutually repulsive force that acts between 
neighboring worlds, which we show in Sec. IV to be 
responsible for a number of generic quantum effects. Note 
also that U N is invariant under translation of the configu¬ 
rations. This leads directly to an analog of the quantum 
Ehrenfest theorem, also shown in Sec. IV. 

While of a simple form, the interworld potential U N in 
Eq. (23) is seen to be a sum of three-body terms, rather than 
the more usual two-body interactions found in physics. 
Further, the corresponding force, acting on configuration 
x n , may be evaluated as 


r N (x n ;X) = -dU N [X)/dx n 

= ^k + tP0-tf»P0], (24) 

which actually involves five worlds, since g„(X) involves 
four worlds: 


<?n(X) = 


1 


1 


{Xfi X n _ i) _ x n-\ | X n X n X n — i 

1 


*„_i - X, 


1 - x„_ 2 _ 


(25) 


(defining x n — — oo for n < 1 and = oo for n > N). 

It is shown in Appendix A that this force corresponds to 
a particular approximation of the Bohmian force r,(x n ), as 
per Eq. (9). 


P'x{ x n ) 


Px( x n+l) — Px( x n) 
X-n+\ 


( 22 ) 


IV. QUANTUM PHENOMENA AS GENERIC 
MIW EFFECTS 


Hence, VP X /P X is approximated by 

« N[P x (x n+ i) -P X ( x n)\ 

r X\ x n) 

1 1 

x n 1 — X n X n — X n _i 

It follows that one may take the interworld potential 
LJ N (X) in Eq. (19) to have the rather simple form 


Here, we show that the MIW approach is capable of 
reproducing some well-known quantum phenomena, 
including Ehrenfest’s theorem, wave packet spreading, 
barrier tunneling, zero-point energy, and a Heisenberg-like 
uncertainty relation. We primarily use the toy model of 
Sec. Ill for this puipose, although a number of the results, 
such as the Ehrenfest theorem and wave packet spreading, 
are shown to hold for broad classes of MIW potentials. 

A. Ehrenfest theorem: Translation invariance 


U N (X) 


h 2 
8 m 


N 


E 

n= 1 


1 

X n +l ~ 


1 I 2 

x n ~ x n— 1 _ 


For any MIW potential U N in Eq. (18) that is invariant 
(23) under translations of the world configurations, such as the 
toy model in Eq. (23), one has 
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d_ 

dt 


x „ = 


1 

m 


(P )n> 


d . , _ . In particular, the position variance for A one-dimensional 

dt' ? w N ’ worlds at time r may be written via Eq. (12) as 


where (^(x^))^ — A -1 <p(x„,p„). Thus, the quan¬ 

tum Ehrenfest theorem corresponds to the special case 
A —► oo, with the more general result holding for any value 
of A. For example, when V(x) is no more than quadratic 
with respect to the components of x, the centroid of any 
phase space trajectory follows the classical equations of 
motion, irrespective of the number of worlds. 

To prove the first part of Eq. (26), note from Eqs. (16) 
and (18) that 




'1 

A 


N 


X 

n= 1 



1 

A 


X-X-(x) 2 n . 


Differentiating with respect to time then gives, using Eqs. (16) 
and (18) with V(x) = 0 an d noting that (d/ df) (p) ;V = Ofrom 
the Ehrenfest theorem (26), 

V N = — [N~ l X • P - {x) N (p) N \ = —Cov Nj (x, p), 
m m 


and 



'P-P 

m 


n=\ 


dU A 


dx„ 


2 

m 2 


(p) 2 n■ 


as required. To prove the second part, note that one also has 
from Eqs. (16) and (18) 

4<p>„ = -<vi 

n= 1 

Now, translation invariance is the condition 

U N (x l + y,..., x N + y) = U N (x l , ...,x N ), (27) 

and taking the gradient thereof with respect to y at y = 0 
yields Yln= t — 0- The second part of Eq. (26) 

immediately follows. 


Now, differentiating the scaling condition [Eq. (28)] 
with respect to A, and setting 2=1, gives 
Yhn=i x n(dU N /dx n ) — -2U n . Hence, recalling that the 
average energy per world is (E) N = A -1 A v , one obtains 

-4 2 

^ = - {E) n - — (p) 2 N = const. 
m m 

Finally, integration yields the variance at time t to be 

2 1 It 2 

T’at(i) = T^(0) 4 Cov N0 (x, p) H 

m m 

(29) 


(p)i 


2 m 


B. Wave packet spreading: Inverse-square scaling 

It is well known that the position variance of a free one¬ 
dimensional quantum particle of mass m, with initial wave 
function 'I' (l . increases quadratically in time [21]: 


2 1 2 1 2 

Var$ x — Var Wn x+ -Cov$ (x,p) H- 

m m 


(P)l„ 




Here, (E)^ and (p) q, denote the average energy and 
momentum, respectively, and Covq,(x,p) denotes 
the position and momentum covariance, (^>\{xp+ 
px)/ 2|T') — (x)^(p)^. This is often referred to as the 
spreading of the wave packet [21]. 

Here, we show that an equivalent result holds for A one¬ 
dimensional worlds, for any interworld potential U ,.y 
satisfying translation invariance and the inverse-square 
scaling property 

U N (Ax u ...,Ax n ) = A~ 2 U n (x u ...,x n ), (28) 

such as the toy model in Eq. (23). 


which is of precisely the same form as the quantum case 
above. Like the generalized Ehrenfest theorem [Eq. (26)] , this 
result holds for any number of worlds A. 

The spreading of variance per se is due to the repulsive 
interaction between worlds having close configurations 
(see previous section). The above result demonstrates that a 
spreading that is quadratic in time is a simple consequence 
of an inverse-square scaling property for the interworld 
potential. Such a scaling is generic, since the quantum 
potential in Eq. (7), which is replaced by U N in the MIW 
approach, also has this property. Note that it is straightfor¬ 
ward to generalize this result to configuration spaces of 
arbitrary dimension K, with V N replaced by the K x K 
tensor (X-X t ) n -(X) n .(X t ) n . Here, the transpose 
refers to the configuration space index k, while the dot 
product refers to the world index n as previously. 

C. Barrier tunneling: Mutual repulsion 

In quantum mechanics, a wave packet incident on a 
potential barrier will be partially reflected and partially 
transmitted, irrespective of the height of the banner. The 
probabilities of reflection and transmission are dependent 
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on the energy of the wave packet relative to the energy of 
the barrier [21]. In the MIW formulation, the same 
qualitative behavior arises as a generic consequence of 
mutual repulsion between worlds having close configura¬ 
tions. This is investigated here analytically for the toy 
model in Sec. Ill, for the simplest possible case of just two 
worlds, N = 2. 

1. Nonclassical transmission 


H 0 (q.q.p.p) 


p 2 p 2 h 2 1 

2 M + 2p + Sp q 2 ' 


with initial conditions p( 0) = 0 and p( 0) = Mv 0 , where 
M = 2m and p = m/2. Thus, the equations of motion for q 
and q decouple, and may be integrated to give 


q(t) =q Q + v 0 t. 



h 2 1 
8 Jiq 2 


h 2 1 

%p q 2 ' 


(31) 


Consider two one-dimensional worlds, the configura¬ 
tions of which initially approach a potential banner V(x) 
from the same side, with kinetic energies less than the 
height of the barrier (Fig. 1). In the absence of an 
interaction between the worlds, the configurations will 
undergo purely classical motion, and so will be unable to 
penetrate the barrier. However, in the MIW approach, the 
mutual repulsion between worlds will boost the speed of 
the leading world. This boost can be sufficient for this 
world to pass through the barrier region, with the other 
world being reflected, in direct analogy to quantum 
tunneling. 

To demonstrate this analytically, consider the toy model 
in Sec. Ill for N = 2. The Hamiltonian simplifies via 
Eqs. (18) and (23) to 


where the 0 subscript indicates initial values. These are 
valid up until one of the configurations reaches the banner 
region. The second equation may be further integrated to 
give the exact solution 

q(t ) = [<?o + {ht/mq 0 ) 2 ] 1/2 . (32) 

Note that this increasing mutual separation is in agreement 
with Eq. (29), with V = q 2 / 4. and is solely due to the 
inverse-square interaction between the worlds in Eq. (30). 

It follows from Eq. (31) that the maximum possible 
separation speed is q „ — fi/(mq 0 ), which can be arbitrar¬ 
ily closely reached if the particles start sufficiently far to the 
left of the barrier. In this case, the speed of the leading 
particle is well approximated by 




(30) 


where V(x) — 0 outside the barrier and has a maximum 
value Vq > 0 in the barrier region. For simplicity, we 
restrict our presentation to the case of equal initial veloc¬ 
ities v Q in the direction of the barrier, with / m < V 0 
(Fig. 1). Thus, in the classical limit h = 0, penetration of 
the barrier is impossible. 

Defining relative and center-of-mass coordinates by 
q •■= x 2 — X{, q ■= (x 2 + X\ )/2, with conjugate momenta 
p = (p 2 ~ P\)/2 and p = p\ + p 2 , respectively, then 
while the configurations of each world are outside the 
barrier region their evolution is governed by the 
Hamiltonian 



FIG. 1. Two one-dimensional worlds, with equal masses m and 
initial speeds v 0 , approach a potential energy barrier of height V 0 . 
Because of mutual repulsion the leading world can gain sufficient 
kinetic energy to pass the barrier. 


x 2 =Q + -qKv 0 + h/ (2mq 0 ) 

at the time it reaches the banner region. Hence, transmission 
through the banner is always possible provided that the 
corresponding kinetic energy is greater than V 0 , i.e., if 

*0 2 m q '' > ^classical \/lV 0 /m, (33) 

where n c]assical is the initial speed required for a classical 
particle to penetrate the banier. 

Equation (33) clearly demonstrates that one world 
configuration can pass through the banner if the initial 
separation between the worlds q 0 is sufficiently small. In 
terms of energy conservation, such a small separation gives 
rise to a conespondingly large interworld potential energy 
in Eq. (30), which is converted into a kinetic energy 
sufficiently large for barrier transmission. Further, the 
second configuration will suffer a corresponding loss of 
kinetic energy, leading to its reflection by the barrier. 

It is seen that even the simplest case of just two 
interacting worlds provides a toy model for the phenome¬ 
non of quantum tunneling, with an interpretation in terms 
of the energy of repulsion between close configurations. 
While this case can be treated analytically, the interactions 
between N > 2 worlds are more complex. It would, 
therefore, be of interest to investigate the general case 
numerically, including tunneling delay times. 
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2. Nonclassical reflection 

The same toy model also captures the nonclassical 
phenomenon that a barrier can reflect a portion of a 
quantum wave packet, even when the incident wave packet 
has a large average kinetic energy [21]. In particular, 
consider the case where the initial kinetic energies of both 
worlds are larger than the barrier height for the “toy” 
Hamiltonian in Eq. (30), i.e., mv\/2 > V 0 . Hence, in the 
classical limit h = 0 both configurations will pass through 
the barrier region. 

It follows from the above analysis that the leading 
configuration is always transmitted—the kinetic energy 
of this world is only further increased by the mutual 
repulsion energy. However, choosing the initial distance 
from the barrier to be large enough for the separation speed 
to approach q 00 , it follows that the second configuration 
will not be transmitted if 


positive. For example, as shown in Appendix B, the 
average energy per world for the toy model in Sec. Ill 
satisfies 

< £ >„><V>„ + (—) jspr. P6) 

This is clearly very similar to the quantum bound in 
Eq. (35), with the latter being precisely recovered in the 
limit N —► oo. 

Remarkably, both the quantum and the MIW lower 
bounds are saturated by the ground state of a one¬ 
dimensional oscillator. The quantum case is well known 
[21]. For the toy model, the corresponding ground-state 
energy is 

(E) N , ground = ( ! - jjj (37) 


V() — 2mq ^ V^o/ m ^classical- (34) 

Indeed, it will not even reach the barrier region if the left- 
hand side is less than zero. This result demonstrates the 
converse to Eq. (33): it is always possible for one 
configuration to be reflected, if the initial separation 
between the worlds is sufficiently small, as a consequence 
of being slowed by the mutual repulsion between worlds. 

D. Zero-point energy: Mutual exclusion 

The classical ground-state energy of a confined system 
corresponds to zero momentum and a configuration that 
minimizes the classical potential V(x). However, the 
corresponding quantum ground-state energy is always 
greater, with the difference referred to as the quantum 
zero-point energy. 

For example, for a one-dimensional quantum system, 
the Heisenberg uncertainty relation (Ax)^(Ap)^ > h/2 
implies 


as demonstrated in Appendix C. Note that this vanishes in 
the classical limit h = 0 (or, alternatively, N = 1), and 
approaches the quantum ground-state energy \ha> in the 
limit N —► oo, as expected. 

E. Heisenberg-type uncertainty relation 

In Sec. IID we noted that the interworld potential in 
Eq. (19) has the form of a sum of nonclassical kinetic 
energies, where the corresponding “nonclassical momen¬ 
tum” of the nth world is given by P" c : = gv(x„;X), with 
the components of g v explicitly defined in Eq. (20). For the 
toy model in Eq. (23), this nonclassical momentum has the 
explicit form 


Pn = 


_Vi+l 


1 -Vi—1. 


(38) 


(defining x 0 = — oo and x N+i = oo, as always). 

It follows immediately that the average of the non¬ 
classical momentum vanishes; i.e., 


(£}*>(V)* + 


h 2 

8/nVa%x 


(35) 


for the average energy of any state 'I' (noting that 
( p 2 )m > Var,|,/o). Hence, for any state with finite average 
energy, the system cannot be confined to a single point, as 
this would imply Va%x = 0. In particular, it cannot be 
confined to the classical ground-state position x min , 
corresponding to the classical ground-state energy 

^min = V(Xmin), alld therefore > ( E )* > V m,n- 

In the MIW formulation, a zero-point energy similarly 
arises because no two worlds can be confined to the same 
position—and to the classical ground-state configuration 
x mm , in particular. This mutual exclusion of configurations 
is a consequence of the repulsion between worlds, which 
forces the interworld potential U N in Eq. (19) to be strictly 


1 N 

V n= 1 

Hence, the energy bound for U N in Eq. (Bl) can be 
rewritten in the rather suggestive form 

(Ax) N (Ap»%> (l-£)^. (39) 

reminiscent of the Heisenberg uncertainty relation for the 
position and momentum of a quantum system. Indeed, a 
nonclassical momentum observable /)',]) may also be 
defined for quantum systems, which satisfies the uncer¬ 
tainty relation A^xAg,/;>^ c > h/2 , and which implies the 
usual Heisenberg uncertainty relation [22]. Thus, Eq. (39) 
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has a corresponding quantum analog in the limit N —> oo. It 
would be of interest to further investigate the role of the 
nonclassical momentum in the MIW approach. 

V. SIMULATING QUANTUM GROUND STATES 

The MIW equations of motion [Eqs. (16) and (18)] in the 
Hamiltonian formulation imply that the configurations of 
the worlds are stationary if and only if 

P« = 0, V x JV(x„) + U N (X)] — 0, (40) 

for all n. In particular, the forces acting internally in any 
world are balanced by the forces due to the configurations 
of the other worlds. 

Unlike quantum systems, the number of stationary states 
in the MIW approach (with N finite) is typically finite. For 
example, for two one-dimensional worlds as per Eq. (30), 
with a classical potential V(x) symmetric about x = 0, 
Eq. (40) implies that the stationary configurations are given 
by x l — —a and x 2 = a, where a is any positive solution of 
V'{a) = h 2 /(\6ma 3 ). For the particular case of a harmonic 
oscillator potential, V(x) — ( I /2)morx 1 . there is just one 
stationary state, corresponding to a — j( h/mco j 1 / 2 . More 
generally, the number of stationary state configurations will 
increase with the number of worlds N. 

The MIW formulation suggests a new approach for 
numerically approximating the ground-state wave function 
and corresponding ground-state energy for a given quan¬ 
tum system. In particular, the ground-state probability 
density is approximately determined by finding the global 
minimum of ^( X J + U (V (X), for suitably large N. 
Below, we give and test a “dynamical” algorithm for doing 
so. First, however, for the purposes of benchmarking this 
algorithm, we calculate the exact ground-state configura¬ 
tions for the toy model of Sec. Ill, for a harmonic oscillator 
potential. 

A. Oscillator ground states: Exact MIW calculation 

For the one-dimensional toy model defined in Sec. Ill, 
with harmonic oscillator potential V(x) = (1/2 )morx 2 , it 
is convenient to define the dimensionless configuration 
coordinates: 

4 — \j2mw/hx n . (41) 

As shown in Appendix C. the unique ground-state con¬ 
figuration for N worlds is then determined by the recur¬ 
rence relation 

( 42 ) 

subject to the constraints 


m 



FIG. 2. Oscillator ground state for N = 11 worlds. The steps of 
the stepped blue curve occur at the values q = ...,fn, 

corresponding to the stationary world configurations 
Xi,x 2 , ...,Xu via Eq. (41). The height of the step between 
and l, +l is P N {£ n ) ■= W '(£„+, - g n )~\ which from Eq. (21) is 
expected to approximate the quantum ground-state distribution 
for a one-dimensional oscillator, P 9 (<*) = for 

large N. The latter distribution is plotted for comparison (smooth 
magenta curve). All quantities are dimensionless. 

^ + ••• + ^ = 0, (£ 1 ) 2 #^- + (6 v ) 2 = iV-l. (43) 

These equations are straightforward to solve for any 
number of worlds. 

As discussed in Appendix C, for small values of N, the 
recurrence relation can be solved analytically, and, in 
general, it can be efficiently solved numerically for any 
given number of worlds N. For example, for N = 11, one 
obtains the configuration depicted in Fig. 2. The corre¬ 
sponding Gaussian probability density for the exact quan¬ 
tum ground state is also plotted as the smooth curve in 
Fig. 2. It is seen that the ground-state configuration of just 
11 worlds provides a good approximation to the quantum 
ground state for many purposes. Note also that the 
corresponding mean ground-state energy follows from 
Eq. (37) as (5/ll)hco « 0.45hco, which is reasonably close 
to the quantum value of 0.5fico. The approximation 
improves with increasing N. 

B. Oscillator ground states: Dynamical MIW algorithm 

The net force acting on each world is zero for a stationary 
configuration of worlds, as per Eq. (40). We exploit this 
fact in the following algorithm to compute a stationary 
configuration. 

Given an arbitrary configuration of N worlds X(0), we 
iterate the following two-step algorithm. (1) Set the 
velocities X(0) to zero, (2) integrate Eq. (11) over a small 
time interval [0, Af], and replace the initial configuration 
X(0) by X(Af). Note that each iteration tends to reduce the 
total energy from its initial value, as the worlds move away 
from maxima of the total potential energy, and the 
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FIG. 3. Convergence of the ensemble of world-particles x n = 
(xi,..., x N ) for N = 11 towards the ground-state configuration as 
a function of the step number in the iteration algorithm described 
in the text. Here, dimensionless positions are used, 
= a JlmwJhx,,, as per Eq. (41), and the temporal step size 
is At = 5 x 10 _2 a> _1 . As the plot shows, convergence is complete 
by step number 6000, at which the distribution of worlds is close 
to the Gaussian quantum ground state. 

velocities are reset to zero in each iteration. The fixed 
points of this iterative map are clearly the stationary states 
of the MIW equations. Hence, after a sufficient number of 
iterations, the configuration will converge to a stationary 
configuration. An alternative algorithm could, instead of 
setting velocities to zero, add a viscous term to the 
evolution, in analogy to the dBB-trajectory-based algo¬ 
rithm of Maddox and Bittner [23]. 

We test this dynamical algorithm for the case of the one¬ 
dimensional harmonic oscillator. We slightly modify the 
interworld force in Eq. (24) by placing auxiliary worlds in 
fixed positions on either side of the configurations 
x { < x x < ■ ■ ■ < x N , with two auxiliary worlds to the far 
left and two to the far right, rather than at — oo and oo. The 
auxiliary worlds have only a tiny effect on the computa¬ 
tions, but ensure the interworld force is well defined for 
n = 1,2..., A. We take advantage of the symmetry of the 
oscillator potential V (x ) to choose an initial configuration 
symmetric about x = 0. 


(E ) N /ti(0 



FIG. 4. Average energy (solid blue line) for the ensemble of 
worlds, as a function of iteration step. Details as in Fig. 3. Note 
that convergence to the exact value for MIW is (1 — l/N)ha>/2 
(red dotted line), which coincides with the quantum mechanical 
value of hi w/2 (green dashed line) for N -> oo. 


Convergence of the algorithm was found to be extremely 
rapid; see Figs. 3 and 4. For example, on a laptop computer 
running a simple Mathematica implementation of the 
algorithm, it takes only 30 sec to converge to the 
ground-state configuration for the case of N — 11 worlds, 
with a corresponding ground-state energy accurate to one 
part in 10 10 in comparison to the exact ground-state energy 
(5/11 )hco following from Eq. (37). 

This dynamical algorithm can be employed to compute 
ground states for all Hamiltonians for which the ground- 
state wave function has a constant phase. The numerical 
advantage over contemporary techniques to compute 
ground-state wave functions is evident when considering 
higher-dimensional problems, i.e., K = DJ > 1. Instead of 
solving the partial differential equation (PDE) H\lf g = 
Eg$> g that lives on a X-dimensional configuration space 
(which memory wise starts to become unfeasible already for 
three particles in three dimensions), one only needs to solve 
for the stationary state of K x N coupled ordinary differ¬ 
ential equations. The quantum distribution T, ; | 2 and 
energy E g of the ground state can then be approximated 
from the stationary distribution of worlds. 

VI. SIMULATING QUANTUM EVOLUTION 

Given the apparent success of the MIW approach as a 
tool for simulating stationary quantum states (Sec. V), it is 
of considerable interest to also investigate whether this 
approach can provide a similarly useful controlled approxi¬ 
mation of the Schrodinger time evolution and, in particular, 
whether it is capable of describing quantum interference, 
which is one of the most striking quantum phenomena. 
Here, we offer evidence that it can, at least for the canonical 
“double-slit” problem. 

In its simplest form, the double-slit scenario comprises the 
free evolution of a one-dimensional wave function for an 
initial value l I' 0 given by a symmetric superposition of two 
identical separated wave packets. We choose real, Gaussian 

6 

4 

2 

q o 

-2 

-4 

-6 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

t 

FIG. 5. A density plot of |4 , ,| 2 for the exact quantum evolution 
of two identical initial Gaussian wave packets in units of the 
initial spread (i.e., a = 1); time is given in units of h/(2m). 
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FIG. 6. Trajectory and (smoothed) density plots of the com¬ 
puted MIW (top) and dBB (bottom) world trajectories for two 
identical initial Gaussian wave packets in units of the initial 
spread (i.e., a = 1) and N = 41 worlds; time is given in units 
of h/ (2 in). 





wave packets with spread a = 1 and initial separation of 4 
units; see Fig 5. As described in Sec. IIB. the corresponding 
initial values for the MIW approach are N worlds, 
X(0) = [jti(0), ...,■%(())], distributed according to |VP 0 1 2 
with zero initial velocities. In order to numerically integrate 
the MIW equations of motion [Eq. (11)] to find the 
corresponding world-particle trajectories X(t), we employ 
a standard Verlet integration scheme [24]. In our example, 
we use N — 41. The computed trajectories X( t) are shown in 
the top two plots of Figs. 6 and 7 using a linearly smoothed 
density and histogram plot. For the comparison with the 
exact quantum solution 'I',, we compute the corresponding 
Bohmian trajectory xj* BB (f) according to Eq. (3) for each 
initial world configuration x, for i— 1,2These 
trajectories are shown in two bottom plots again using a 
linearly smoothed density and histogram plot. Note that, 
due to equivariance, the Bohmian configurations xf iR (t) are 
distributed according to | 2 for all times t as the xf BB (0) 
are distributed according to 'I' 0 | 2 . Therefore, the bottom 
density plot in Fig. 6 is a discretized plot of |'I',| 2 , as shown 
in Fig. 5 (this discretization, corresponding to using just 
N = 41 Bohmian worlds and linear smoothing, is respon¬ 
sible for the blurriness relative to Fig. 5). 


FIG. 7. Smoothed histogram plots of the computed MIW (top) 
and dBB (bottom) world trajectories for two identical initial 
Gaussian wave packets in units of the initial spread (i.e., a = 1) 
and N = 41 worlds; time is given in units of h/(2m). 

By comparison of the plots, we conclude that the MIW 
approach is at least qualitatively able to reproduce quantum 
interference phenomena. Furthermore, it should be stressed 
that with respect to usual numerical PDE techniques on 
grids the MIW approach could have several substantial 
computational advantages: As in dBB, the worlds X{t) 
should remain approximately |T' r | 2 distributed for all times 
t, which implies that without increasing N, the density is 
naturally well sampled in regions of high density even 
when the support of 'I' ( becomes very large. The computa¬ 
tional effort is then automatically focused on high-density 
regions of configuration space, i.e., on the crucial regions 
where numerical errors have to be minimized to ensure 
convergence in the physically relevant L 2 norm. 
Furthermore, transmission and reflection coefficients 
can be computed simply by counting worlds. These 
advantages also hold true for methods of integration in 
the Schrodinger equation on comoving grids as, e.g., 
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proposed by Wyatt [20]. However, there are reasons to 
expect that the MIW approach may be numerically more 
stable (see Ref. [25] for a discussion of such numerical 
instabilities in conventional approaches) and less perfor¬ 
mance intensive. 

There is much work to be done on the question of time 
evolution using our MIW approach. Even in the example 
above, the convergence as a function of N, the number of 
worlds, remains to be investigated. Also, there is the issue 
of wave functions with nodes (or even finite regions of 
zero values). These allow for different wave functions 
T'o(g) having identical Po(q), and identical VS Q (q) 
everywhere that / J 0 (r/) ^ 0. From Sec. IIB, these different 
wave functions will all generate the same initial 
conditions for our ensemble of many worlds, for any 
finite value of N, and thus are indistinguishable in the 
MIW approach, contrary to quantum predictions (see 
also Sec. X of Ref. [12]). For example, the wave 
functions < T > o(<z) and |rf> 0 (ry) | are of this type, where 
<f> 0 (g) is any real wave function that takes both positive 
and negative values. 

However, it is plausible that restricting the moments of 
the quantum energy to be finite would typically eliminate 
all, or all but one, of such wave functions, by enforcing the 
continuity of V'T 0 everywhere [26,27], Finally, beyond the 
toy model, there is the question of whether the general 
approach of many interacting worlds allows successful 
numerical treatments of more general situations, with 
common potentials and K = DJ > 1. 

VII. CONCLUSIONS 

We introduce an approach to quantum phenomena in 
which all quantum effects are due to interactions between a 
large but finite number N of worlds, and probabilities arise 
from assigning an equal weighting to each world. A 
number of generic quantum effects, including wave packet 
spreading, tunneling, zero-point energies, and interference, 
are shown to be consequences of the mutual repulsion 
between worlds (Sec. IV). This alternative realistic inter¬ 
pretation of quantum phenomena is also of interest in not 
requiring the concept of a wave function. The wave 
function does not exist other than as an epiphenomenon 
in the notional limit that the initial distribution of worlds 
approaches |T' 0 (x)| 2 and the initial velocity of each world 
approaches Eq. (3) as N -> oo. 

For finite N, our MIW approach can only ever give an 
approximation to quantum mechanics, but since quantum 
mechanics is such an accurate theory for our observations, 
we require this approximation to be very good. In this 
context, it is worth revisiting the question raised in 
Sec. IIC: What restrictions must one place on the initial 
distribution of world configurations and velocities so that 
nearly quantum behavior (that is, observations consistent 
with quantum mechanics for some wave function) will be 
experienced by macroscopic observers in almost all worlds? 


The answer may depend on the details of the interworld 
potential, but we can suggest directions to explore. 

As discussed in Sec. IIC, it may be necessary to impose 
the quantum velocity condition of Eqs. (3) and (4) using 
some smoothly varying (on the scale of the separation 
between nearby worlds) single-valued complex function 
on configuration space. However, it is conceivable that 
other (perhaps even most, by some measure of typicality) 
initial conditions would relax, under our interacting-world 
dynamics, to conditions approximating quantum behavior, at 
least at the scale that can be probed by a macroscopic 
observer. This is one of the big questions that remains to be 
investigated. One might think that this idea would never 
work, because the velocities would relax to some type of 
random Maxwellian distribution. However, this may well not 
be the case, because the many-body interaction potentials 
and forces in our MIW approach—which generically repro¬ 
duce various quantum phenomena as per Secs. IV-VI —are 
very different from those assumed in classical statistical 
mechanics (see also the “gas” example discussed in Sec. I A). 

Other matters for future investigation include how spin 
and entanglement phenomena such as teleportation and 
Bell-inequality violation are modeled in the MIW approach. 
The latter will require studying the case of worlds with 
configuration spaces of at least two dimensions (correspond¬ 
ing to two one-dimensional systems), and will also allow 
analysis of the quantum measurement problem (where one 
system acts as an “apparatus” for the other). This may help 
clarify the ontology and epistemology of any fundamental 
new theory based upon the MIW approach to quantum 
mechanics. 

In the context of entanglement, it is worth comparing our 
MIW approach with conventional many-worlds approaches. 
The latter are often motivated by the desire, first, to restrict 
reality to only the wave function, and second, to avoid the 
explicit nonlocality that arises from entanglement in other 
realist versions of quantum mechanics. Our approach is, by 
contrast, motivated by the desire to eliminate the wave 
function. It furthermore elevates the nonlocality of quantum 
mechanics to a kind of “supemonlocality”: particles in 
different worlds are nonlocally connected through the 
proposed MIW interaction, thus leading, indirectly, to 
nonlocal interactions between particles in the same world. 

Turning from questions of foundations and interpreta¬ 
tions to applied science, the MIW approach provides a 
promising controlled approximation for simulating quan¬ 
tum ground states and the time-dependent Schrodinger 
equation, as discussed in Secs. V and VI. In particular, we 
show that it is able to reproduce quantum interference 
phenomena, at least qualitatively. Quantitative comparisons 
with different initial conditions, convergence as a function 
of the number of worlds N, and generalizations to higher 
dimensions, are a matter for immediate future work. 
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APPENDIX A: CONNECTION OF BOHMIAN 
FORCE WITH THE HAMILTONIAN 
FORMULATION OF MIW 

We briefly note here how the interworld potential U N in 
Sec. IIC is related to a corresponding approximation of the 
quantum potential Q in Sec. II A. We also exhibit a direct 
correspondence between the interworld force r w (q; X f ) 
and the Bohmian force r r (q) for the case of the toy ID 
model defined in Sec. III. 

For any particular wave function 'I', (q), the dBB 
equation of motion [Eq. (5)] for the world-particle is 
generated by the Hamiltonian 

p) ■■=^2 { -f^+V(x) + Q t (x), 

k= l zm 


following from integration by parts and Eq. (7) [2]. In 
particular, since the N world-particles have configurations 
sampled according to | X P(q) | 2 , it immediately follows that 


1 

N 


72=1 


i N K 

v££ 


n= 1 k= 1 


n 2 

%m k P 2 


ppy 

dq 1 ) 


9k= 


Hence, U N in Eq. (Al) may alternatively be replaced by the 
form in Eq. (19), corresponding to a suitable approximation 
of P and its derivatives in terms of the configurations of the 
N worlds. Here “suitable” means that we must check that 
the corresponding MIW force is given by r, v (x„; X) = 
—V x JJ N (X) in the limit N —► oo. 

For the toy MIW model of Sec. Ill, we can demonstrate 
directly that the force in Eq. (24) is an approximation of the 
Bohmian force, as per Eq. (9) of Sec. II. To do so, we first 
note that, for the one-dimensional case, the derivative of 
any regular function q>(x) can be approximated at x = x n 
by either of 


(p'(x n ) 


<p(x„+i) - <p(x n ) ~ q>(x„) - cp(x n _ i) 
Xn+ 1 — %n X n — X n — 1 


for sufficiently large N. Hence, using Eq. (21), one obtains 


where <2(q) is the quantum potential corresponding to 
\D,(q) [identifying the momentum component p k with the 
right-hand side of Eq. (3)]. Hence, the evolution of N 
Bohmian world-particles, each guided by H/^q), is pre¬ 
cisely described via the time-dependent Hamiltonian 


tf*'(X,P) P») 


n= 1 

NK/„k\2 N 


_ \ ^ (Pn 

‘4-' f—f 2 m k 

n= 1 k= 1 




72=1 


72=1 


Thus, if Q t ( q) (and its gradient) can be approximated 
by some function Q(q;X f ), depending on the configura¬ 
tions X, of the N worlds [assumed to sample |'I',(q)| 2 ], 
then one immediately has a suitable corresponding 
Hamiltonian formulation of the MIW approach, with an 
interworld potential U N in Eq. (18) of the alternative 
form 


1 d' 

P{x n ) dx 


cp(x n )xND + cp(x n ) 


ND_cp(x n ), 


(A2) 


where the difference operators D ± are defined by 


D + q> n := q> n+l - (p n , D_cp n ■■= cp n - cp n _ x 


for any sequence {((>,, }■ 

Now, as may be checked by explicit calculation from the 
definition of the quantum potential in Eq. (7), for a one¬ 
dimensional particle, the Bohmian force in Eq. (6) can be 
written in the form 

r, = (ti 2 /Am) (1 / P) [P(P'/P)']' - , (A3) 


i.e., 


r M) = 


h 2 

4m 


1 d' 

Pin) 2 

1 d' 

P{q) dq 


P{q)dq 


p (q). 


N 

f/ w (X) = ^e(x„;X) 

72=1 


(Al) 


where the derivative operators act on all terms to their right. 
Applying Eq. (A2) then gives 


to that in Eq. (19). 

To show how the above form is related to the positive- 
definite form in Eq. (19), consider the relation 

/<**l*(q)l 2 e(q) = (f £) 2 


r t (x n )K—N 2 D + lP(x n ) 2 D + D_P(x n )} 


h 2 

4m 4 


1 




- 2 d + d_ 


x„ - x, 


n- 1. 


= r N {x n \X), 


(A4) 
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as desired, where the second line uses Eq. (21), and the 
last line follows by expansion and direct comparison 
with Eq. (24). 


APPENDIX B: LOWER BOUND FOR THE 
INTER WORLD POTENTIAL 

To obtain the bound in Eq. (36), first let /i, ...,/# and 
g 0 ,g l ,g N be two sequences of real numbers such that 
9o — 9 n = 0- It follows that 


N 

(E) n = N~ 1 H n > N~ l Y V(x n ) + N-'UnW 

n= 1 


> -mco 2 V N + 


N - 1 


N- 1 fico 


N 

2NmwV N 


2 h 2 1 
8 m V N 
(N-l)h 


N 


4 


(N-l)h 


N - Ihco 

> -. 

~ N 2 


2NmcoV N 


N 

Yfn(9n 9n—\) 

n =1 


N -1 

TM n +1 “ fn)9n > 

n =1 


as may be checked by explicit expansion. Further, the 
Schwarz inequality gives 


r N 


L«=l 


r N 


nr N 


^ 'j fn(.9n 9n— l) — ^ ) (/';/) ~ ^ ^ (9 n 9n— l) 


Li—1 


-I Ll=l 


with equality if and only if f n = a(g n — p„_i) for some a. 
Combining these results then yields 


The second line follows via (x 2 ) N > V, v and inequality 
(Bl), and the last line via z + 1/z > 2. 

To show that this inequality chain can be saturated, as per 
Eq. (37), the conditions for equality must be checked at 
each step. It follows that the lower bound is achievable if 
and only if (i) the momentum of each world vanishes, i.e., 
Pi = ■ ■ ■ = Pn = 0, and (ii) the configuration satisfies 
both 


{x) N = 0, 


V N — 


N - 1 h 
N 2mm 


(Cl) 


and the Schwarz inequality saturation condition 




n= 1 


~9n- 1) 2 > 


[E^(/» + i-/»)gJ 

E Llifn ) 2 


a 


x„ = 


x n +1 X n 


x„ -X, 


n— 1 


(C2) 


The particular choices f n —x n - (x) N and g n = 1/(x„ +1 —x„) 
(with x 0 = —oo and x N+ \ = oo) immediately yield, using 
the definition of U N in Eq. (23). the lower bound 

„ , rf KSiP (w-i) 2 tf l 
" _ 8 m E?.i(*«-W ») 2 N s 1 1 

for the interworld potential energy, with equality if and only 
if the Schwarz inequality saturation condition 


from Eq. (B2), for some constant a (using (x) N = 0). These 
conditions generate a second-order recurrence relation with 
fixed boundary conditions, yielding a unique solution for 
any given number of worlds N. The corresponding unique 
ground-state configurations converge to the quantum 
Gaussian ground-state wave function for N oo, as 
discussed in Secs. IIC and VA. 

To determine the ground-state configuration, for a given 
number of worlds N, note that summing each side of 
Eq. (C2) from n = 1 up to any n < N gives 


x„ - x „ = 


%n +1 x n 


X„ — X 


(B2) 


n -1 


is satisfied. Equation (36) then follows via Eq. (18), as 
desired. 

Similar bounds can be obtained for other choices of the 
interworld potential, but are not discussed here. 


APPENDIX C: GROUND-STATE ENERGY 
AND CONFIGURATION OF TOY 
MODEL OSCILLATOR 

To derive the ground-state energy in Eq. (37), for the ID 
toy model oscillator, note first that for a harmonic oscillator 
potential V(x) = (1 /2)ino> 2 x 2 , one has 


1 X] + • • • + x„ 


x n +1 x n 


a 


(C3) 


This is also satisfied for n = N (defining x^+i = oo 
as usual), since x, + • • • + x N = 0 from Eq. (Cl). 
Further, noting that x x < ■ ■ ■ < x N and (xj + • ■ ■ + x n )/n < 
(x) = 0, it follows from Eq. (C3) that a < 0. Hence, since 
summing the squares of each side of Eq. (C2) over n gives 


8 m a 2 

NV n = -^U n , 


and inequality (Bl) is saturated by the ground state, it 
follows via Eq. (Cl) that 
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NV n _ _ h_ 

N — 1 2 mm 


(C4) 


Defining the dimensionless configuration coordinates 
~ (— a)~ x ! 2 3 4 x n = \/lmm/fix n , (C5) 

conditions (Cl), (C2), and (C3) then simplify to 

£i 4-§* : &v = 0 , {£\Y - 1 , (C6) 

fn+l=£n~ . , 1 , . • (C7) 

Cl -r-T Cn 


Since these equations are invariant under —► — g n , the 

symmetry property 


Zn — ^N+l—n (C8) 

follows from the uniqueness of the ground state. 

For N G {2,3,4, 5}, the recurrence relation [Eq. (C7)] 
reduces to solving an equation no more than quadratic in 
(£i Y- allowing the ground state to be obtained analytically. 
For example, for N = 3, one finds 

& = -l, & = 0, ft = l. 

while for N = 4, one obtains 

e _ V / 7+Tl7_ e 

Sl ~ 

£2 = li + 2 = -£3. 

More generally, it may be used to express ■ ■ ■. ‘S[,v/ 2 i in 
terms of £ 1; where [z] denotes the integer part of z. The 
ground-state configuration may then be numerically deter¬ 
mined by solving the condition (fi) 2 + • • • + (£[;v/ 2 ]) 2 = 
(N — l)/2 for <*! [following from Eq. (C6) and the above 
symmetry property]. 
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